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Abstract. The generation and evolution of magnetic fields in the plasma 
accreting into a rotating black hole is studied in the 3+1 split of the Kerr 
metric. Attention is focused on effects of the gravitomagnetic potential. 
The gravitomagnetic force appears as battery term in the generalized 
Ohm's law. The gravitomagnetic battery is likely to saturate at higher 
field strength than the classical Biermann battery. 

The coupling of the gravitomagnetic potential with electric fields 
appears as gravitomagnetic current in Maxwell's equations. In the mag- 
netohydrodynamic induction equation, this current re-appears as source 
term for the poloidal magnetic field, which can produce closed magnetic 
structures around an accreting black hole. In principle, even self-excited 
axisymmetric dynamo action is possible, which means that Cowling's 
anti-dynamo theorem does not hold in the Kerr metric. 

Finally, simulations of the aO. dynamo in accretion flows into the 
hole are presented. I assume a simple expression of a in this relativistic 
context. The modes of the dynamo are oscillating for dynamo numbers 
which are typical for accretion disks. In a zero angular momentum flow 
into a Kerr black hole there is still shear, i.e. the angular velocity Vt of the 
plasma equals the angular velocity of space, uj, and it has been speculated 
(Meier 1998) that even then a dynamo could operate. This is shown to 
be unlikely due to the rapid accretion. 



Introduction 



The influence of the gravitomagnetic field of a rotating compact object on elec- 
tromagnetic fields has been studied for some 25 years (Wald 1974, Ruffini & 
Wilson 1975, Blandford &: Znajek 1977). The coupling of the gravitomagnetic 
potential with a magnetic field results in an electromotive force. Currents driven 
by this electromotive force may extract rotational energy from a black hole. This 
energy could power relativistic jets by Poynting flux (but see also, e.g., Punsly 
1996). Cast in the language of the 3+1 split of the Kerr metric. Maxwell's equa- 
tions, together with the 'ingoing wave boundary condition' for electromagnetic 
fields at the horizon, led to the Membrane Paradigm (Thorne et al. 1986). Since 
a black hole does not carry its own magnetic field, not to mention kCauss fields 
required for the Blandford-Znajek process to be efficient, strong magnetic fields 
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must either be accreted into the black hole from the outer accretion disk, or 
have to be generated and amplified in the plasma surrounding the black hole. 

The generation of magnetic fields by a battery operating in the plasma close 
to a rotating black hole was studied by Khanna (1998b). It was shown that the 
gravitomagnetic force may play a crucial role in the battery. Khanna & Camen- 
zind (1996a) studied the possibility of an axisymmetric gravitomagnetic dynamo 
(or uj^l dynamo), in which the coupling between the gravitomagnetic potential 
and an electric field is a source for the poloidal magnetic field. This theoretical 
result invalidates Cowling's anti-dynamo theorem (see also Niihez 1997), but 
self-excited growing dynamo modes could not yet be numerically varified for 
simple kinematics (Khanna & Camenzind 1996b). Egi et al. (1998) have given 
a criterion for growing modes of the dynamo, i.e. that the Poynting flux 
from close to the horizon (extracted rotational energy of the hole) be positive. 

Recently, Meier (1998) speculated that, in a zero angular momentum flow 
into a Kerr black hole, alternatively to the cjfi dynamo, an acu dynamo might 
operate. My simulations of such a sceanrio show that, due to the relativistic 
accretion velocity, an auj dynamo is not very likely. 

Section ^ gives an introduction to the derivation of MHD in the 3-1-1 split 
of the Kerr metric. The generalized Ohm's law for an electron-ion plasma is 



presented in Sec. 2.1.. The gravitomgnatic battery (along with the relativistic 
equivalent of Biermann's battery) are discussed in Sec. 2.2\ In Sec. ^ the MHD 
induction equation in the 3-1-1 split of the Kerr metric is presented. Applications 
are the gravitomagnetic dynamo (Sec. 3.1. and |3.2.| ) and, in Sec. |3.3.| , the ail 



dynamo, or the auj dynamo, respectively. Throut the paper I set G = 1 = c. 



2. The MHD description of an electron-ion plasma 

The formulation of MHD requires the relativistic definition of a plasma as center- 
of-mass fluid of its components (Khanna 1998a). The plasma is assumed to be 
a perfect fluid and is defined by the sum of the ion and electron stress-energy 
tensors, which contain a collisional coupling term: 

x=i,e 

Subscripts i, e refer to ion and electron quantities, respectively. Superscripts de- 
note the rest-frame in which the quantity is defined, where refers to the plasma 
rest-frame. In the 3-1-1 split (into hypersurfaces of constant Boyer-Lindquist time 
t, filled with stationary zero angular momentum fiducial observers) T°^^ splits 
into the total density of mass-energy e and momentum density S 

e = (Pm +P'^'^)7^ ~ /5m7^ S = {p'm+p')l'^v^ p'^-i'^v (2) 

and the stress-energy tensor of 3-space with metric h 

T= {p^+ p')j'^v 1^ V + p' p^^'^v igi V + p' h . (3) 

The approximate expressions hold for a 'cold' plasma. Charge density and cur- 
rent density are given by 

Pc = Pci + Pee = Zeriiji - eUeJe j = ji + je = ZeUiJiVi - eUeJeVe ■ (4) 
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All quantities resulting from the split are measured locally by FIDOs. 



2.1. The generalized Ohm's law in the 3+1 split of the Kerr metric 

In the 'cold' plasma limit, the local laws of momentum conservation for each 
species can be re-written as equations of motion, which can then be combined 
to yield the generalized Ohm's law for an electron-ion plasma 

j Zriiyi ^ j^B V(agP°) 47r7e Pc'iv 
— « E + —vxB-- + ^ ^ + ^T^Pc9 + ^-^-^ 

(5) 
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with the conductivity a = e^n^jra^Vf. = Wpg/47rfc as measured in the plasma 
rest frame, the electron plasma frequency Wpe, the factor of gravitational red- 
shift ag (with g = — Vlnog) and the gravitomagnetic tensor field H= a~^V/3 . 

P = P^e,/, = —oje^ is the gravitomagnetic potential, which drags space into differ- 
ential rotation with angular velocity uj. Note that, in the single fluid description, 
the gravitomagnetic force drives currents, only if the plasma is charged in its rest 
frame. Tp is the proper time in the plasma rest frame. The derivation requires 
the assumption that the species are coupled sufficiently strong that their bulk 
accelerations 



dry, 



+ 7x U'x • V 



Og dt 



(7x^x) (6) 



are synchronized. The same is required for the gravitomagnetic accelerations, 

i.e. I H ■{'j^Vi)- H ■{'^lve)\ < I H -(7^ If the MHD-assumption of "synchro- 
nized accelerations" is not made. Ohm's law contains further current acceleration 
terms, inertial terms and gravitomagnetic terms (Khanna 1998a), which may be 
important for collisionless reconnection and particle acceleration along magnetic 
fields. This topic will be discussed elsewhere. 

In the limit of quasi-neutral plasma {Zni ~ rie) and 7c ~ 71 ~ 7, Eq. (^) 
reduces to 

J ^ a^{E + vxB)-—{jxB) + -^V{agpl) , (7) 

which contains all the terms, familiar from the non-relativistic generalized Ohm's 
law, but no gravitomagnetic terms. 

2.2. The gravitomagnetic battery 

The generation of magnetic fields by a plasma battery was originally devised by 
Biermann (1950) for stars. He showed that, if the centrifugal force acting on a 
rotating plasma does not possess a potential, the charge separation owing to the 
electron partial pressure cannot be balanced by an electrostatic field, and thus 
currents must flow and a magnetic field is generated. 
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In Khanna (1998b) I have re-formulated Biermann's theory in 3+1 spht 
of the Kerr metric. The base of this battery theory is Ohm's law of eq. (|^). 
Assuming that electrons and ions have non-relativistic bulk velocities in the 
plasma rest frame, superscripts i, e, can be dropped. With p = pi + = 
{rii + ne)kT, the impressed electric field (lEF), E^'^^ = V {agPe) / enjag, can be 
re-expressed with the aid of the equation of motion for a 'cold' quasi-neutral 
plasma to yield 

f^m rrii / ^ t; . ^^ d(jv)\ Z (j x B + {v-j)E 



(y^U K^9+H ijv)-^^]+ ^ (8) 

(Z + 1)6 V ar / [Z + l)enj 

T is the proper time in a FIDO frame; i.e. d/dr^ = 'jd/dr. The criterium for 
magnetic field generation is that V x agE^^^ ^ . Here I restrict the discussion 
to the gravitomagnetic lEF Eg^ . The function part of OgE'gm is 

(/3 • V + V/3 •) {^v) = [p\-iv% + ^if\H^ e,- 

= -7?;<^cD2Vcu-a;((7t,'-),<^e, + (7?;<^),^eV) , (9) 

where a) = (hif,^)^^'^ and | denotes the covariant derivative in 3-space. In ax- 

isymmetry OgEg^ is clearly rotational, unless some freak 7 should manage to 
make 'yv'^oj'^ a function of cu alone. Thus the gravitomagnetic force drives a 
poloidal current and generates a toroidal magnetic field. Only if v'^ is non- 
axisymmetric, the gravitomagnetic lEF drives toroidal currents. The total IFF 

(cKg^gm + ckg-^ciass) Hkely to rotational in general. This will be quantified for 
specific velocity fields elsewhere. 

In presence of a weak poloidal magnetic field the Biermann battery is limited 
due to modifications of the rotation law by the Lorentz force, rather than by 
ohmic dissipation. Then the contribution of the centrifugal force to the IFF 
becomes irrotational already at weak toroidal fields (Mestel & Roxburgh 1962). 
The gravitomagnetic battery term, on the other hand, is only linearly dependent 
on V. The equilibrium field strength should therefore be higher than for the 
Biermann battery. 



3. The MHD induction equation in the 3+1 split of the Kerr metric 

In this section I review the axisymmetric dynamo equations in the 3+1 split of 
the Kerr metric (Khanna & Camenzind 1996a). Ohm's law is assumed to be of 
the standard form for a quasi-neutral plasma; Hall-term and IFF are neglected. 
Combining Maxwell's equations (Thorne et al. 1986) with Ohm's law yields the 
MHD induction equation 

^ = V X (^{agv X fi) - ^ (V X (agB) + (^p • Vu;)we^)^ + (^p • Va;)cDe^ . 

The term standing with the magnetic diffusivity 77 is the current density, which, 
via Ampere's law, contains the coupling of the gravitomagnetic field with the 
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electric field. In axisymmetry this is simply the shear of the poloidal electric 
field in the differential rotation of space, to. Another induction term is the shear 
of the poloidal magnetic field by uj. This generates toroidal magnetic field out 
of poloidal magnetic field even in a zero-angular-momentum flow. 

3.1. The gravitomagnetic dynamo 

Introducing the flux ^ of the poloidal magnetic field and the poloidal current T 
^ii = ^ J Bp-dA = OjA^ ^ " ^ / "s-^P ■ " OgwS'^ , (11) 
where A'^' is the toroidal component of the vector potential, eq. (|To| ) splits into 
^4- + a,(.-p . V)^ - ^.^V.. . W) - . f ^ 



jag 



Tvp - ^VT ) X e: 



Vuj (12) 



^ + ag(t;p-V)r + agi;2;^(v.||)-ag(^2^.(^Vr) 

= agcD(V^' X e^) • VVL . (13) 

These equations are the relativistic equivalent of the classical axisymmetric dy- 
namo equations. It is important to note, however, that no mean-field approach 
was made, but ^ has source terms anyway. They result from ■ Vlo. Obvi- 
ously, Cowling's anti-dynamo theorem does not hold close to a rotating black 
hole. Growing modes of this gravitomagnetic dynamo were shown to exist for 
steep gradients of the plasma angular velocity O (Nuiiez 1997). According to Egi 
et al. (1998) growing modes require that the Poynting flux carrying rotational 
energy extracted from the hole be positive. For simple accretion scenarios, grow- 
ing modes could not be found in kinematic numerical simulations (Khanna &: 
Camenzind 1996b). If, however, magnetic field is replenished by an outer bound- 
ary condition, the gravitomagnetic source terms generate closed loops around 
the black hole. 

3.2. The magnetic field structure in the accretion disk close to the 
hole 

In an accretion disk, magnetic fields may be advected into the near-horizon area, 
where gravitomagnetic effects may become important. This can be simulated 
by advection/diffusion boundary conditions for F = ^ , T 

1^ + = ^ + ^F.., , (14) 

on r] on rj 

where d/dn is the derivative along the outer boundary normal. Figure |^ shows 
the stationary final state of a time-dependent simulation (with turbulent mag- 
netic diffusivity) , in which |-Bp^out ,out \ — 1/50. For such a dominantly toro- 
idal magnetic field the gravitomagnetic source terms are strong enough to change 
the topology of ^. This may influence the efficiency of the electromagnetic ex- 
traction of rotational energy from the hole. 
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Figure 1. Left: Contours of magnetic flux ^ showing a quadrupo- 
lar magnetosphere of an accreting, rapidly rotating black hole (a = 
0.998M, 

Qvisc — 0.25). Right: The gravitomagnetic current as source 
of ^. The disk is marked by long-dashed lines. Solid contours corre- 
spond to positive values, short-dashed contours indicate negative val- 
ues. The range of contours is given below the boxes. Time is measured 
in diffusive timescales (see below), = GM/c^. 



3.3. The aft dynamo in the Kerr metric 

In the innermost region of an accretion disk around a black hole there may also be 
a turbulent source term of a-type. Without any knowledge of the physical source 
of the term (convection or magnetic shear instability) or its mathematical form 
in the relativistic context, one can try to assess the physical regime (magnetic 
diffusivity, accretion velocity, rotation law) in which growing modes of an afl 
dynamo exist. For a simple mean-field ansatz (Khanna & Camenzind 1996a) 
the equations of the kinematic aQ dynamo are identical to Eqs. ([l^ ) and ([l^ ) 
augmented by the a-source term, aT, for the flux and r/ replaced by r/turb- In 
analogy to the expression for a in classical disks, I assume 

a = {a^Rolln/H) f{z)/f{H/2) = (Saga^iscHQ) f{z)/f{H/2) , (15) 

where H is the disk scale height, Ro is the Rossby number and Ovisc is the vis- 
cosity parameter of standard accretion disk theory. The factor ag is added in 
order to suppress the source close to the horizon, where the accretion velocity ap- 
proaches the speed of light (in properly derived mean-field equations there would 
probably be a "Rossby number" correlated to the accretion velocity instead). 
The vertical dependence of a is modelled with f{z) = tanh{z) eyip[—{z/H)'^]. 
The turbulent diffusivity is described as 

??turb = avisc-f^^^^w/r sin(6') , (16) 

with a vertical scaling (exp[— (z/i7)^] -|- 0.1)/1.1. The boundary condition at 
r = 10 Tg is d'^ /dn = and T = 0. Figure |2| shows the first part of a simulation 
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Figure 2. A simulation of an ail, dynamo with symmetric initial cur- 
rent in a quasi-Keplerian accretion disk. Kerr parameter a = 0.998M, 
"vise = 0.065. Solid contours have positive values, dashed contours 
have negative values. Kinks at the outer edge are artefacts of trans- 
forming data from the spherical grid into a carthesian plot. Simulation 
continued in Fig. |3|. 
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Figure 3. The simulation of Fig. ^ has reached a slowly growing, 
oscillating eigenmode with a period of ~ 3tdifr- Deviations from equa- 
torial (anti-)symmetry are probably due to insufficient resolution in 
^-direction. Not shown here: The inclusion of non-linear a-quenching 
leads to severe symmetry breaking and chaotic behavior. 
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R/r^ R/r^ 

-0.0151452953, -0.001 / 0.001, 0.02079677396 -0.01222794596, -0.001 / 0.001, 0.004189541098 



Figure 4. Simulation with same setup as above, except that the ac- 
cretion velocity at r ^ 2 rg is 10 times higher. 

with an initial current T, which is symmetric with respect to the equatorial plane, 
and \I' = 0. The parameters are the Kerr parameter a = 0.998 M, avisc = 0.065, 
and the angular momentum of the accreting plasma is 99.999% of the Keplerian 
value at r > r^s and constant within, which yields an accretion velocity of 
^ 0.003c at r = 3rg, increasing to c at the horizon. The dynamo is in a slowly 
growing quadrupolar mode, oscillating with a period of about three diffusive 

timescales tdiff = r^/r/o ~ 2 10^ sec Mg {^)~^ {^)~^ ■ 

The same setup, but with lower angular momentum in the accretion disk 
(99.9% of the Keplerian value), which corresponds to a radial velocity of ~ 0.03c 
at r = 3rg, is in a decaying mode, which demonstrates that accretion impedes 
dynamo action (Fig. |^. 

3.4. auj dynamo action in a zero-angular momentum flow? 

It was mentioned above that the shear of space does also induce a toroidal 
magnetic field (cf. eq. fl^). In eq. ( |l3| ) this shear term is obscure, but still 
there, hidden in (V^' x e^) • oc Bp ■ Vfi = Bp ■ 'V^Ogv'^ + uj). In a zero- 
angular-momentum flow v"^ = (or, equivalently Q = lo) and thus the current T 
is solely generated by the shear of space. Moreover, Vw is significantly steeper 
than Vr^K, which means that Bp ■ Voj is a strong source term. 

Meier (1998) speculated that, alternatively to the gravitomagnetic dynamo 
described above, there could be an aoj dynamo in a zero-angular momentum 
accretion flow, with a being due to the magnetic shearing instability. Such a flow, 
however, accretes at relativistic velocities (;^ 0.1c at r = 10 rg), which should 
suppress any dynamo action. This conclusion is supported by the simulation 
shown in Fig. |5|. The ^' and T loops in the corona are transient and depend 
on the description of -q and a (here described as in Eq. |]l5|, but not suppressed 
by ag in order to have an upper estimate of the source term). Parameters are 
Ovisc = 0.003 and H{rh) = 0.8 rg. A wider parameter study is in progress. 
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t = 20.078000 ^ t = 20.078000 T 




-0.4619527161, -0.001 / 0.001, 0.08680785447 -1.476681352, -0.001 / 0.001, 1.476729631 



Figure 5. Simulation of an auj dynamo. The dynamo generates tran- 
sient structures in the corona close to the horizon. Within the disk 
there are no signs whatsoever of dynamo action. The field and current 
are completely determined by the relativistic advection. T is shown in 
logarithmic contours. 
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